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RESUMEN 

El resumen sera traducido al espanol por los editores. There are several lines of evidence which suggest 
that the relativistic outflows in gamma-ray bursts (GRBs) are collimated into narrow jets. The jet structure 
has important implications for the true energy release and the event rate of GRBs, and can constrain the 
mechanism responsible for the acceleration and collimation of the jet. Nevertheless, the jet structure and its 
dynamics as it sweeps up the external medium and decelerates, are not well understood. In this review I discuss 
our current understanding of GRB jets, stressing their structure and dynamics. 

ABSTRACT 

There are several lines of evidence which suggest that the relativistic outflows in gamma-ray bursts (GRBs) are 
collimated into narrow jets. The jet structure has important implications for the true energy release and the 
event rate of GRBs, and can constrain the mechanism responsible for the acceleration and collimation of the 
jet. Nevertheless, the jet structure and its dynamics as it sweeps up the external medium and decelerates, are 
not well understood. In this review I discuss our current understanding of GRB jets, stressing their structure 
and dynamics. 
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1. INTRODUCTION 

Gamma-ray bursts (GRBs) are produced by a 
highly relativistic outflow from a co mpact sourc e 
(for a comprehensive recent review see IPiranll2005h . 
Early GRB models featured a spherical outflow, 
mainly for simplicity. However, other astrophysical 
sources of relativistic outflows such as active galactic 
nuclei and micro-quasars are in the f orm of narrow 
bipolar jets. One might argue (e.g., lRhoadslll997h 
that, in analogy to other such sources, GRBs might 
also be collimated into narrow jets. 

The initial Lorentz factor during the prompt 
gamma-ray emission is very high, Tq > 100, and 
therefore we observe emission mainly from very small 
angles, < T^ 1 < 10 -2 rad, relative to our line of 
sight. This is a result of relativistic beaming (i.e. 
aberration of light), an effect of special relativity, 
which causes an emission that is roughly isotropic 
in the rest frame of the emitting fluid (as is gener- 
ally expected under most circumstances) to be con- 
centrated mostly within an angle of around its 
direction of motion in the lab frame, where r » 1 
is the Lorentz factor of the emitting fluid in the lab 
frame. For this reason, the prompt gamma-ray emis- 
sion probes a region of solid angle ~ ttTq 2 , or a frac- 
tion - Tq 2 / 4 ~ iO -7 - 10~ 4 - 5 of the total solid an- 
gle, and cannot tell us whether the outflow occupies 
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a larger solid angle. 

Therefore, more direct evidence in favor of jets 
in GRBs had to await t he discovery of af terglow 
emi ssion in the X-ray llCosta et all Il997t) . opti- 
cal llyan Paradiis et al.lll997|) . and radio l|Frail et alJ 
11997^ 7 that lasts for days, weeks, and months, respec- 
tively, after the GRB. The afterglow is believed to 
be synchrotron emission from the shocked external 
medium. As the relativistic outflow expands out- 
wards it sweeps up the surrounding medium and 
drives a strong relativistic shock into it, called the 
forward shock, while the ejecta are decelerated by 
a reverse shock. Eventually, most of the energy 
is transfcred to the shocked external medium be- 
hind the forward shock, and t he flow approaches a 
spher ical self-similar evolution felandford fc McKeel 
Il976j) , gradually decelerating as it sweeps up the ex- 
ternal medium. 

This is valid not only for an initially spherical 
outflow, but also for the interior of a jet with an 
angular size larger tha n T^ 1 (as appears to be th e 
case for GRB jets, e.g. IPanaitescu fc Kumarll2002(l . 
as long as r -1 remains smaller than the angular size 
of the jet and the interior of the jet is out of causal 
contact with its edges (i.e. before the jet break time) . 
Therefore, before the jet break time the isotropic 
equivalent energy of the jet, E{ so , is relevant (both for 
its dynamics and for the resulting emission), while 
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at very late times as the jet becomes sub-relativistic 
and approaches spherical symmetry its true energy, 
E, is relevant. In the intermediate regime things are 
more complicated, and are discussed in this review. 

The forward shock is responsible for the long 
lived afterglow emission, while the reverse shock pro- 
duces a shorter lived emission, that peaks in the 
optical or NIR on a time-scale of tens of seconds, 
when the reverse sh ock crosses the shell of ejecta (the 
"optical flash". e.g.lAkerlof et ajll999USa,ri fc Piranl 
Il999albl : iMeszaros fc R.eeslll999|) . The shocked out- 
flow gradually cools adiabatically and the peak of 
its emission shifts to lower frequencies, until after 
about a day it peaks in the radio (the "radio flare" 
Kulkarni et aJ]ll999btlFrail et allEooil iBerger et all 



2003a|) . During the afterglow, the Lorentz factor T 



of the emitting shocked external medium decreases 
with time as it accumulates more mass, causing the 
visible region of 9 < T _1 around the line of sight 
to increase with time. This enables us to probe the 
structure of the outflow over increasingly larger an- 
gular scales. 

Different lines of evidence suggest that the rel- 
ativistic outflows in GRBs arc collimated into nar- 
row jets. A compelling, although somewhat indi- 
rect, argument comes from the very high values for 
the energy output in gamma rays assuming isotropic 
emission, E 1> - ISOi that are inferred for GRBs with 
known redshifts, z, which approach and in one case 
(GRB 991023) even exceed M Q c 2 . Such extreme en- 
ergies in an ultra-relativistic outflow arc hard to pro- 
duce in models involving stellar mass progenitors. If 
the outflow is collimated into a narrow jet (or bipolar 
jets) that occupies a small fraction, /j, < 1, of the 
total solid angle, then the strong relativistic beaming 
due to the very high initial Lorentz factor (ro > 100) 
causes the emitted gamma rays to be similarly col- 
limated. This reduces the true energy output in 
gamma rays by a factor of to E 1 = fbE 7 ^ so , 
thus significantly reducing the energy requirements. 

Estimates of the energy in the afterglow 
shock from late time radio observations when 
the flow is only mildly relativistic and starts to 
approach spheric al symmetry (often called "ra- 
dio calorimetry" : iFrail. Waxman fc Kulkarnil 2000t 
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ically yield Ek ~ 10 515 erg which lends some sup- 
port for the true energy being significantly smaller 
than -E 7i iso- One should keep in mind, however, that 
these are only approximate lower limits on the true 
afterglow energy, and the latter can in principle be 
much higher (see, e.g.. lEichler fc WaxmanlEoOfil) . 
Furthermore, there is good (spectroscopic) ev- 



idence that at least some GRBs of the long-soft 
class occur together (to within a fe w days) with 
a core collapse supernov a of Type Ic l|Stanek et all 
l2003UHiorth et al.ll2003h . In such cases the average 
Lorentz factor must be (r) < 2 for a spherical explo- 
sion, since the accreted mass does not significantly 
exceed the ejected mass, and only a fraction of the 
rest energy of the former can provide the kinetic en- 
ergy for the latter. Therefore, only a small fraction of 
the ejected mass can reach T > 100 which is required 
in or der to power the GRB, and hydrodynamic anal- 
ysis llTan. Matzner fc McKeel l200lt IPerna fc Vietrl 
120021) shows that it would carry a small fraction of 
the total energy which is insufficient to account for 
the high end of the observed values of E 7 ^ so . For a 
jet the ejected mass can be much smaller than the 
accreted mass so that (r) ^> 1 is possible, in addi- 
tion to the smaller E-y that is implied by the same 
observed i? 7 ,i S0 . 

A more direct line of evidence in favor of 
a narrowly collimated outflow comes from achro- 
matic breaks seen in the afterglow light curves of 



1999aHHarrison et al!l999HStanek et al Jl999ll200ll 


Bcrger et all 12000 


: lHalDern et al.ll2000l IPrice et all 


200ll 


ISagar et al.l 


200lllJensen et alJl200lf>. In fact. 



such a "jet break" in the afterg low light curve was 
predicted before it was detec ted (jRb.oad£Hl997Lll999l 
ISari. Piran fc Halpernlll"999(l . The cause of the jet 
break in the light curve is discussed in M2.6I 

The properties of GRB jets are of fundamental 
importance since they pertain to the GRB energy re- 
lease, event rate, and the progenitor model through 
its ability to produce a particular jet structure. In 
particular, a good understanding of the jet struc- 
ture and dynamics are crucial in order to reliably 
address these vital issues. This review focuses on 
the dynamics of the jet as it sweeps up the exter- 
nal medium and decelerates O and on its angular 
structure (fl, stressing the constraints that may be 
derived from various observations. The conclusions 
are discussed in [Jlj 

2. THE JET DYNAMICS 

This section begins by presenting three different 
approaches to the calculation of the jet dynamics, in 
order of increasing complexity: semi-analytic models 
(§ 12.1(1 . simplifying the dynamical equations by inte- 
grating over the radial profile of the jet (S I2.2|) . and 
full hydrodynamic simulations f§ l2.3fl . The main re- 
sults of the different approaches are described and 
compared. Next (§ 12. 4[) there is a brief description 
of the typical assumptions that are made in order to 
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calculate the afterglow emission. The afterglow im- 
age is discussed in § 12. 51 along with potential methods 
for resolving it and constraining its angular size, as 
well as how its morphology and the evolution of its 
size may help us learn about the jet dynamics and 
the external density profile. Finally, the cause of the 
jet break in the afterglow light curve is discussed in 

2.1. Simple Semi- Analytic Models 

The first approach that had been adopted 
for calculating the jet dy namics was using a 
simple semi- analytic model (|Rhoadsl Il997l Il999j) . 
Many different variat ions on this basic approach 
have followed (e.g.. ISari. Piran fc Halpernl Il999t 
IPanaitescu fc_ Meszarosl Il999t iKumar fc Panaitescvj 

isT l20ori 



Moderski. Sik ora fc Buli 



^UU; uyioclerski^jikQ] 

lOren. Nakar fc Piranl |2004 



For simplicity we 
present h ere an a nalysi s that largely follows the 
model of iRhoadsl l)l999|) . and which captures the 
main features of this type of models. 

The basic underlying model assumptions are (i) 
a uniform jet within a finite half-opening angle 9j 
with an initial value 9o that has sharp edges, (ii) 
the shock front is part of a sphere at any given lab 
frame time and the emitting fluid behind the forward 
shock has a negligible width, (iii) the outer edge of 
the jet is expanding sideways at a velocity c s ~ c in 
the local rest frame of the jet, (iv) the jet velocity 
is always in the radial direction and 9j <C 1. Under 
these assumptions, the jet dynamics are obtained by 
solving the ID ordinary differential equations for the 
conservation of energy and particle number. 2 The 
lateral expansion velocity in the comoving frame, c s , 
is usually identified with the sound speed, in which 
case c s » c/\^3 while the jet is relativistic. However, 
this does not have to be the case: it could in principle 
be either much smaller (c s <C c), or as large as the 
thermal speed (i.e. c s sa c wh ile the jet is relativistic; 
ISari. Piran fc Halpernlll999T) . 

The lateral size of the jet, R±, and its radius, i?, 
are related by R± ss 9jR. We have 



dR±_ w 6jdR + c s dt' 



(i) 



2 For the adiabatic energy conserving evolution considered 
here, the equation for momentum conservation is trivial in 
spherical geometry, and does not constrain the dynamics. For 
a narrow (9j <g; 1) highly relativistic (r 2> 1) jet, the equation 
for the conservation of linear momentum in the direction of the 
jet symmetry axis is almost identical to the energy conserva- 
tion equation. When the jet becomes sub-relativistic the con- 
servation of energy and linear momentum force it to approach 
spherical symmetry, and once it becomes quasi-spherical then 
again the momentum conservation equation becomes irrele- 
vant. 



where dt' = dt^/T sa dR/cT and 

dR~ R\dR V ~ cRT{R) 



(2) 



Eq. |21 suggests that 9j ~ 9 + c s /cT, and there- 
fore the jet expands significantly when T drops to 
~ c s /c8q. This can occur after the edge of the 
j et becomes visible (when T ~ for c s < c 

l|Panaitescu fc Meszaroslll999|) . Once the jet begins 
to expand sideways significantly, then to zeroth order 
6, j cx r _1 and therefore energy conservation suggests 
that R — const, since E — r 2 9 2 R 3 p 0Xt (R)c 2 . Here 
p cxt = AR~ k is the external density, which is as- 
sumed to be a power law in radius. 3 As is shown 
below, a more careful analysis shows that T9j slowly 
decreases with radius (Eq. I10JI while 6j grows very 
rapidly with radius (Eq. I13|) . 

The total swept-up (rest) mass, M(R), is accu- 
mulated as 

H a 2^i?) 2 Pc *t(i?) - 2irAR 2 - k 6*(R) , (3) 

where the factor of 2 is since a double sided jet is 
assumed. As long as the jet is relativistic, energy 
conservation takes the form E f=a r 2 Mc 2 , which im- 
plies that Md{T 2 ) = -T 2 dM, and 



dR 



r dM 

2M~dR 



= -nAR 2 - k 9 2 (R) 



r(fl) 

M(R) 



(4) 



One can numerically integrate equations J5J, and 
© thus obtaining 9j(R), M(R), and T{R). Alterna- 
tively, one can use the relation E ps T 2 Mc 2 (energy 
conservation) which reduces the number of free vari- 
able to two, and solve equations (J2J and Q. Chang- 
ing variables to a dimensionless radius, R = R/Rj, 
where 



gives 



d6i_ 
dR 
dT 
dR 



Rj ~ I nAc 2 



0s 



l/(3-fc) 



RY{R) 



— 2 r>2 — k n2 



(R)T 3 (R) 



(5) 

(6) 
(7) 



The initial conditions at some small radius Rq <C 1 
(just after the deceleration radius) are 9j(Ro) = 6q 
and 



T(Ro) 



- k P s 5 -(3-fe)/2 

~9- R ° 



(8) 



3 We consider here and throughout this review only k < 3 
for which the shock Lorentz factor decreases with radius for a 
spherical adi abatic blast wave during th e self-similar stage of 
its evolution iBlandford fc McKedll973) . 
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Fig. 1. The jet dynamics according to the simple semi- 
analytic model that is described in the text (solid lines), 
that is obtained by numerically solving equations © and 
JTj with the initial conditions given by Oj(Ro) = do and 
equation JSJ. We have used j3 a = c a /c = 3 _1//2 , which 
corresponds the the sound speed of a relativistically hot 
fluid, and show results for a uniform external medium 
(k = 0) and for a stellar wind (k = 2). Also shown are 
the analytic approximations for Ra cc < R < Rj (dashed- 
dotted lines) and for Rj < R < Rnr (dashed lines, ac- 
cording to equations 1101 . 1131 and 1141 with b = 1/4). 
For T9j at Rj < R < -Rnr we also show (by the dotted 
line) the higher order approximation given in footnote [I] 



Equations © and J7J imply, 



dk 



R 



R 



2-fc 



-(r%) 3 



(9) 



If one assumes that the first term becomes dominant 
at R > 1 then this equation implies T9j ss (3 S In R, 
which in turn implies that the second term would be 
dominant, rendering the original assumption incon- 
sistent. The same applies if the opposite assumption 
is made, that the second term is dominant (in this 
case T9j oc i?( fe_3 )/2 which implies that the first term 
would be dominant). This implies that the two terms 
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Fig. 2. The jet Lorentz factor T as a function of the 
observed arrival time of photons emitted along the line 
of sight Ti os « J dR/2cT 2 , for the simple semi-analytic 
model illustrated in Figure (solid lines). Both T and 
ri os are normalized to their values at Rj extrapolated 
from R <C Rj (Tj and 2I OSj j, respectively). Also shown 
are the asymptotic scalings at T os <C Tj and Ti os 3> Tj . 



must remain comparable, and therefore 4 

TOj m /3 s R- {3 - k)/3 ■ 



(10) 



A similar conclusion can be reach by taking the ratio 
of equations JBJ and Q which implies that 

d(T- 3 )=j3- 3 R 3 - k d{9 3 j ) . (11) 

Substituting equation (|l(Jfl into equation © 
yields 



dR 



and 



?o cxp 



.(3-fc) 

r » A^-3)/3 exp 



9 3 R~ k/3 



R(3-k)/3 
3 



(3-fc) 



_R(3-*0/3 



(12) 



(13) 



(14) 



where & « 1/4 is determined numerically. 

The results of this simple semi-analytic model are 
illustrated in Figures Q and [3 In practice, the dy- 
namical range between the onset of the exponential 
lateral spreading of the jet (at Rj) and the non- 
relativistic transition (at i?NR.) is quite limited. This 

4 In order to satisfy equation J^J another term with a 
smaller power in R is required, F6j a f3 a {R( k ~ s ^ 3 _|_ 
jj2(fc-3)/3^3 _ k)/9], but only the leading term is shown in 
equation 1101 , T his result i s consistent with equation (8) of 
iKumar fc Panaitescul (2000) where the second term in that 
equation dominates in the relevant regime. 
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fact is ignored in these figures, and a wide dynamical 
range is shown in order to better isolate the charac- 
teristics of this intermediate stage (Rj < R < Rnr)- 
The dynamical transition at R ~ Rj is much more 
gradual for a wind environment (fc = 2) compared 
to a uniform density medium (fc = 0). This leads 
to a much smoother and m ore gradual jet break 
in th e afterglow light curve (|Kumar fc Panaitescul 
which would be hard to detect. 
The results d erived here ar e somewhat differ- 
ent than those of iRhoadsl (^)99) who obtained V oc 
cxp(— R) and 9j oc R^ 1 exp(i?) at Rj < R < -Rnr 
for a uniform external medium (fc = 0), an d they 
are clo ser (though not identical 5 ) to those of IPiranl 
(|2000|h This demonstrates the sensitivity of such 
semi-analytic models to the exact assumptions that 
are made. Nevertheless, despite the differences in 
their details, all of these semi-analytic models for 
the jet dynamics share a similar main prediction: a 
very fast lateral expansion (where the jet half open- 
ing angle 9j typically grows exponentially with the 
radius R) after the jet break time. As is discussed in 
M2.2I and W2.'S\ more detailed numerical calculations 
of the jet dynamics, which better capture the rele- 
vant physics, contradict this result and show that 
the degree of lateral expansion is very modest as 
long as the jet is relativistic. Therefore, a simple 
and useful approximation for (semi-) analytic calcu- 
lations would be that the jet does not expand side- 
ways altogether, retaining its original opening angle 
and evolving as if it were part of a spherical flow, as 
long as it is relativistic. 

2.2. Intermediate Approach: Integrating over the 
Radial Profile 

The over-simplified treatment of the jet dynam- 
ics in simple semi-analytic models, and the fact that 
different such models obtained different results, put 
into question the validity of those results and moti- 
vated more careful studies of the jet dynamics. A 
proper treatment of this problem requires a full hy- 
drodynamic simulation (in at least 2D) and is dis- 
cussed in the next subsection. However, since such 
simulations arc very challenging numerically, an in- 
termediate approach between simple semi-analytic 
models and full hydrodynamic si mulations can be 
useful . This was attempted by iKumar fc Granotl 
(|200ril) and is briefly described here. Under the as- 
sumption of axial symmetry, the dynamical equa- 
tions are reduced to two spatial dimensions. The 

5 The difference arises since there it was assumed that 
M(R) oc p elct (R)R 2 L R, while here the differential form is used, 
dM oc p c ^ t {R)R\dR. 




Fig. 3. The dynamics of a jet with an initially (at 
a lab frame time to) Gaussian profile where the en- 
ergy per solid angle e and the Lorentz factor Y are 
given by e, T — 1 oc exp(— 8 2 /28 2 ), calculated using a 
scheme where the dynamical equations are simplified by 
integrating over the rad ial profile of the shocked fluid 
iKumar fc Granotjliooll . The parameters used in this 
calculation are C = 0.035 rad, F(6 = 0,t ) = 200, 
e(0 = 0,t ) = 10 53 /4tt erg/sr, n ex t = 10 cm" 3 . The 
different panels show the evolution of V — 1, e, and the 
comoving lateral velocity v' g /c = Tvg /c. The different 
curves are for different lab frame times: v' g /c increases 
with time, while e and V — 1 decrease with time at the 
center of the jet and increase with time at the sides of 
the jet behind what appears as a shock in the lateral 
direction which propagates to larger angles 0. 

dynamical equations can be reduced to a set of 
one dimensional (ID) partial differential equations 
(PDEs), and thus greatly reducing the computa- 
tional difficulty involved in solving these equations, 
by integrating over the radial profile of the shocked 
fluid in the jet. The motivation for this, apart from 
making the calculations much easier, is that most of 
the shocked fluid is concentrated within a very thin 
layer behind the shock, of width A ~ i?/I0T 2 in the 
lab frame (i.e. the frame of the unperturbed exter- 
nal medium). This suggests that integrating over the 
radial profile might not introduce a very large error 
in the calculation of the jet dynamics. 

The results of this method are shown in Figure 
for a jet with an initial Gaussian profile in T — 1 
and in the energy per solid angles e = dE/dQ. Such 
a jet can be thought of as a smoother and there- 
fore more realistic version of a uniform jet, since it 



6 



GRANOT 




1 10 

e 

Fig. 4. The dynamics of a "structured" jet where initially 
r-1 = 200/(1 + e 2 /6 2 c ) and e = e /(l - 6 2 /6 2 c ) with e = 
10 53 /4tt erg/sr, 6> c = 0.02 rad, and n oxt = 10 cm" 3 . The 
format is similar to Figure |3] Again v' e jc increases with 
time while V — 1 decreases with time. For a structured jet, 
as opposed to an initially Gaussian jet, a shock does not 
develop in the lateral direction and e(9) remains almost 
unchanged as long as the jet is relativistic. 

has a roughly uniform core and relatively sharp (but 
still smooth) wings. A shock appears to develop in 
the lateral direction, because of the very steep initial 
angular profile in the wings of the Gaussian. Never- 
theless, the lateral expansion remains modest as long 
as the jet is relativistic. This can be seen both from 
the small (compared to c) lateral velocity in the co- 
moving frame, and from the fact that e(9) does not 
change very much compared to its initial distribu- 
tion. The modest degree of lateral spreading is in 
stark contrast with the results of semi-analytic mod- 
els. 

Figure 0] shows the resulting dynamics for a 
"structured" jet (which is discussed in S|3J| where 
r and e are initially power laws with the angle 6 
from the jet symmetry axis, outside of some narrow 
core. Again, there is very little lateral expansion (i.e. 
the comoving lateral velocity remains <C c and e(9) 
hardly deviates from its initial profile) as long as the 
jet core is relativistic. 

2.3. Hydrodynamic Simulations 

The most reliable method for calculating the jet 
dynamics is using hydrodynamic simulations. This 
is a formidable numerical task for the following rea- 
sons. First, it requires a hydrodynamic code with 




Fig. 5. A 3D view of a relativ istic impulsive jet at the 
last time step of the simulation dGranot et all200iTl . The 
outer surface represents the shock front while the two in- 
ner faces show the proper number density (lower face) 
and proper synchrotron emissivity (upper face) in a log- 
arithmic color scale. 

special relativity that is accurate over a large range 
in the four velocity u — Tf3, from it k T > 1 to 
u k, (3 <§C 1. Second, the shocked fluid in the jet is 
concentrated in a very thin layer behind the shock, 
of width A ~ i?/10r 2 in the lab frame, which is 
extremely narrow at early times when r > 1, and 
therefore very hard to resolve properly. 

More specifically, from considerations of causal- 
ity, significant lateral expansion could in princi- 
pal occur when T becomes comparable to 0q , so 
that ideally one would want to start with an ini- 
tial Lorentz factor To 3> Qq 1 , an d in practice 
we need at least To^o ^ a few. Observed jet 
break times suggest 0.05 < 9q < 0.2 and there- 
fore require r > 20 - 100. If 100A 2 cells arc 
needed in order to resolve the shell of width A ~ 
i?/10T 2 then the minimal cell size in the initial 
time step (denoted by the subscript '0') needs to 
be of the order of S ~ lO -6 ^" 1 (r /30)- 2 i? ~ 
10- 6 A 2 " 1 (r 6lo/3)- 2 (6lo/0.1) 2 J Ro- The minimal num- 
ber of cells required to resolve the initial shell is 
A min ~ 9 R ( AoS- 2 ~ lO 7 A 2 (r o o /3) 2 (0o/O.l)- 1 . 
The total number of cells in each time step can be 
Atot ~ A m in if the code uses adaptive mesh re- 
finement (AMR). Otherwise, for a fixed cell size, 
A tot /A rain > i? /A ~ lO 4 (r o o /3) 2 (VO.l)- 2 . 

Because of the numerical difficulty involved, very 
few attempts have been m ade so far ijGranot et alJ 
l200ll ICannizzo et al.lE00^) . In the following I shall 
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Fig. 6. The jet radius (upper panel) , Lorentz factor minus 
one (middle panel), and angle 9 from the jet symmetry 
axis (lower panel), as a function of the lab frame time 
tiah ~ R/c (in day s), from a hydrodynamic simulation 
(iGranot et al.ll200ll) . 



concentrate on the results of IGranot et alJ l)200l|) . 6 
The initial conditions were a cone of half-opening 
angle f5 = 0.2 rad, taken out of the spherical 



6 The calculations of ICannizzo et alJ 120041) suffer from 
poor numerical resolution. 



iBlandford & McKed l)l976|) self-similar solution with 
an (isotropic equivalent) energy of i?k,iso = 10 52 erg 
and a uniform external density of n ex t = 1 cm~ 3 . 
The initial Lorentz factor of the fluid just behind the 
shock was Tq rs 16.8 corresponding to TqOq « 3.4. 

The results of the simulation are illustrated in 
Figures El and El While the number density does 
not change significantly between the front and the 
sides of the jet, the emissivity is large only at the 
front of the jet, within its initial half-opening angle 
(6 < 9q), and drops sharply at 8 > 9q. This causes 
the emissivity weighted values of the Lorentz factor 
and radius to be close to their maximal values (which 
are also obtained at the front of the jet). While 
the sides of the jet contribute a small fraction of 
the total emissivity, their emission can dominate the 
observed flux for lines of sight outside the initial half- 
opening angle of the jet, as discussed in H3.5I The 
overall egg-shaped structure of the shock front is very 
different from the quasi-spherical structure assumed 
in ID semi-analytic models. Moreover, the degree 
of lateral expansion is very modest as long as the 
head of the jet is relativistic, in contradiction with 
the very rapid lateral expansion predicted by semi- 
analytic models. 

2.4. The Afterglow Emission 

The dominant emission mechanism during the af- 
terglow stage is believed to be synchrotron emission. 
This is supported by the detection of linear polar- 
ization at the level of ~ 1% — 3% in several opti- 
cal or NIR afterglows (see H3.3fl . and by the shape 
of the broad band spectrum, which consists of sev- 
eral power-law segments that smoothly join at some 
typical break frequencies. Synchrotron self-Compton 
(SSC; the inverse- Compton scattering of the syn- 
chrotron photons by the same population of relativis- 
tic electrons that emits the synchrotron photons) can 
sometimes dominate the afterglow flux i n the X-rays 
l|Sari fc Esinll200li lHarrison et ai]l200ll) . 

It is usually assumed that the electrons are 
(practically instantaneously) shock-accelerated into 
a power law distribution of energies, dN/dje cx j~ p 
for 7 e > j m , and thereafter cool both adiabatically 
and due to radiative losses. Furthermore, it is as- 
sumed that practically all of the electrons take part 
in this acceleration process and form such a non- 
thermal (power-law) distribution, leaving no ther- 
mal com ponent (which is not at all clear or justi- 
fied; e.g. lEichler fc WaxmanlEoO^ . The relativistic 
electrons are assumed to hold a fraction e e of the in- 
ternal energy immediately behind the shock, while 
the magnetic field is assumed to hold a fraction es 
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Fig. 7. The afterglow synchrotron spectrum, calculated 
for the iBlandford fc McKee) lll976f) spherical self-similar 
solution, under standard assumptions, using the accurate 
form of the synchrotron spectral emissivity and integra- 
tion over the emission from the whole volume of shocked 
materia l behind the forward (afterglow) shock (for de- 
tails see iGranot fc Sar3l200^ . The different panels show 
the five possible broad band spectra of the afterglow syn- 
chrotron emission, each corresponding to a different or- 
dering of the spectral break frequencies. Each spectrum 
consists of several power law segments (PLSs; each shown 
with a different color and labeled by a different letter 
A-H) that smoothly join at the break frequencies (num- 
bered 1-11). The broken power law spectrum, which 
consists of the asymptotic PLSs that abruptly join at 
the break frequencies (and is widely used in the litera- 
ture), is shown for comparison. Most PLSs appear in 
more than one of the five different broad band spectra. 
Indicated next to the arrows are the temporal scaling of 
the break frequencies and the flux density at the differ- 
ent PLSs, for a uniform (ISM) and stellar wind (WIND) 
external density profile. 

of the internal energy everywhere in the shocked re- 
gion. This is a convenient parameterization of our 
ignorance regarding the micro-physics of relativistic 
collisionless shocks, which are still not sufficiently 



well understood from first principals. 

The spectral emissivity in the co-moving frame 
of the emitting shocked material is typically approxi- 
mated as a broken power-law (in some cases the more 
accurate fun ctional form of the synchrotron emission 
is use d, e.g. IWiiers fc Galama1ll999t IGranot fc Saril 
120021) . Most calculations of the light curve assume 
emission from an infinitely thin shell, which repre- 
sents the shock front (some integrate over the volume 
of the shocked fluid taking into acco unt the appro- 
priat e radial profile of the flow, e.g. IGranot fc Saril 
I2002L see Figure |7|). One also needs to account for 
the different arrival times of photons to the observer 
from emission at different lab frame times and loca- 
tions relative to the line of sight, as well as the rel- 
evant Lorentz transformations of the emission into 
the observer frame. SSC is included in some (not 
all) works, although it can also effect the synchrotron 
emission through the enhanced radiative cooling of 
the electrons. 

2.5. The afterglow Image 

The apparent surface brightness distribution and 
size evolution of the afterglow image on the plane 
of the sky can potentially provide very useful infor- 
mation about the structure and dynamics of GRB 
jets, as well as about the radial dependence of the 
external density. Furthermore, polarimetry (or even 
spectral polarimetry) of a resolved afterglow image 
could provide valuable information on the magnetic 
field structure behind collisionless relativistic shocks, 
which is not well understood theoretically. However, 
most GRBs are at cosmological distances (z > 1) 
and the angular size of their afterglow image is of 
the order of a micro-arcsecond (/xas) after a day or 
so, making it extremely difficult to resolve the image. 

During the self-similar spherical evolution stage 
(before the jet break time, for a jet), the afterglow 
image has circular symmetry around the line of sight 
(where the surface brightness depends only on the 
distance from the center of the image), and is con- 
fined within a circle on the sky with a radius 

{f \V 8 n \ 5 / 8 
3-91 (ff) (fe) (* = 0) 

/„ \l/4/ f \3/4 
2-39 (If) (^) (ft = 2) 

(15) 

(see Figure [HJ where S52 is the isotropic equivalent 
afterglow kinetic energy in units of 10 52 erg, i days 
is the observed time in days, and the external den- 
sity is assumed to be a power law with the distance 
R from the central source, p C xt = n cxt m p = AR~ k , 
where hq = n oxt /(l cm~ 3 ) for a uniform external 
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observed 
image 



r(R. 
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Fig. 8. Schematic illustration of the equal arrival time 
surface (thick black line) , namely, the surface from where 
the photons emitted at the shock front arrive at the same 
time to the observer (on the far right-hand side). The 
maximal lateral extent of the observed image, R±, is lo- 
cated at an angle , where the shock radius and Lorentz 
factor are R* and T* = r s h(i?»), respectively. The area 
of the image on the plane of the sky is S± = -kR\. The 
shock Lorentz factor, r s h, varies with radius R and angle 
6 from the line of sight along the equal arrival time sur- 
face. The maximal radius Ri on the equal arrival time 
surface is located along the line of sight. If, as expected, 
r s h decreases with R, then Ti = Y s ^{Ri) is the minimal 
shock Lorentz factor on the equal arrival time surface, 
(from iGranot. Ramirez-Ruiz fc Loebll2005l) . 

density (k = 0), and A* = A/{5 X 10 11 g cm" 1 ) for 
a wind-like external density profile (fc = 2) as might 
be expected for a massive star progenitor. This cor- 
responds to an angular radius of 



R± 

d A 



1.61 uas 

d,A,27.7 

0.98 ^as 

d-A,27.7 



1/8 

n J \ 1+z J 



( Em. V /4 ( tday3 

\A.) { 1+z 



5/8 



3/4 



(fc = 0) 



(fc = 2) 



(16) 

where cIa{z) is the angular distance to the source, 
and d,A,27.7 is d,A m units of 10 27 7 cm f=a 5x 10 27 cm. 7 

More generally, the afterglow image size during 
the self-similar spherical stage scales with the ob- 
served time as R± cx £( 5_fc )/ 2 ( 4_fc ) . The image size 
grows super-luminally with an apparent expansion 
velocity of r s h(i?*)c. The expected afterglow images 
in this self-similar regime are shown in Figures El and 
I1UI The normalized surface brightness profile within 
the afterglow image is independent of time due to 
the self-similar dynamics, and changes only between 
the different power law segments of the synchrotron 
spectrum, and for different external density profiles. 
The image becomes increasingly limb-brightened at 
higher frequencies, and for smaller values of k. 

Below the self-absorption frequency the spe- 



7 For a standard cosmology (Qm = 0.27, Qj± = 0.73, h 
0.72) d A (z) has a maximum value of 5.37x 10 27 cm (rfA,27.7 
1.07) for 2 = 1.64. 




Fig. 9. The afterglow images for different power law seg- 
ments of the spectrum, for a uniform ( k = 0) and wind 
(k = 2) external dens ity profile (from IGr anot fc Loebl 
1200 ll . calculated for the lBlandford fc McKedll976l sphe~ 
ical self similar s olution, using the formalism of 
IGranot fc Sarill200^ . Shown is the surface brightness, 
normalized by its average value, as a function of the 
normalized distance from the center of the image, r = 
Rsm8/ R± (where r = at the center and r — 1 at the 
outer edge). The image profile changes considerably be- 
tween different power-law segments of the afterglow spec- 
trum, F v cx v 13 . There is also a strong dependence on the 
density profile of the external medium, p cxt cx R~ k . 



cific intensity (surface brightness) represents the 
Rayleigh-Jeans portion of a black-body spectrum 
with the blue-shifted effective temperature of the 
electrons at the corresponding radius along the front 
side of the equal arrival time surface of photons to 
the observer (i?» < R < Ri in Figure |SJ). Above 
the cooling break frequency the emission originates 
from a very thin layer behind the shock front where 
the electrons whose typical synchrotron frequency is 
close to the observed frequency have not yet had 
enough time to significantly cool due to radiative 
losses. This results in a divergence of the surface 
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Fig. 10. An illustration of the expected afterglow image 
on the plain of the sky, for thre e different power law 
segme nts of the spectrum (from lOranot. Piran fc Saril 
1999allrl. assuming a unif orm external density and the 
Blandford fc McKed 19761) self-similar solution. The im- 
age is more limb brightened at power law segments that 
correspond to higher frequencies. 

brightness at the oute r edge of the image llSa,rilll998l 
IGranot fc Loebll200ll) . 

After the jet break time the afterglow image is 
no longer symmetric around the line of sight to the 
central source for a general viewing angle (which is 
not exactly along the jet symmetry axis), and its 
details depend on the the hydrodynamic evolution 
of the jet (so that in principal it could be used in 
order to constrain the jet dynamics). Therefore, a 
realistic calculation of the afterglow image during the 
more complicated post-jet break stage requires the 
use of hydrodynamic simulation, and still remains to 
be done. 

The afterglow image may be indirectly resolved 
through gravitational lensing by a star in an inter- 
vening galaxy (along, or close to, our line of sight 
to the source). This is since the angular size of the 
Einstein radius (i.e. the region of large magnifica- 
tion around the lensing star) for a typical star at a 
cosmological distance is ~ 1 /zas (hence the name 
micro- lensing) , and therefore comparable to the af- 
terglow image size after a day or so. Since the af- 
terglow image size grows very rapidly with time, dif- 
ferent parts of the image sample the regions of large 
magnification (close to the point of infinite magnifi- 
cation just behind the lensing star) with time, and 
therefore the overall magnification of the afterglow 
flux as a function of time probes the surface bright- 
ness profile of the afterglow image. This results in a 
bump in the afterglow light curve which peaks when 
the limb-brightened outer part of the image sweeps 
past the lensing star, where the peak of the bump 
is shar per the more limb-br ightened the afterglow 
image l)Granot fc LoeblboOlh . It has been suggested 
that an achromatic bump in the afterglow light curve 
of GRB 000301C after ~ 4 days might have been du e 
to micro- lensing l|Garnavich. Loeb fc Stanekll2000() . 
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Fig. 11. Tentative fits to the constraints on the image size 
(or diameter 2R±) of the radio afterglow of GRB 030329 
at different epochs, for different external density profiles 
and different assumptions on the lateral spreading of the 
jet. The physical parameters and external density pro- 
file for each model are indicated (the viewing angle is 
along the jet symmetry axis). Model 1 features relativis- 
tic sideways expansion in the co-moving rest frame of 
the je t material, while model 2 has no late ral spreading 
(from IGranot. Ramirez-Ruiz fc Loebll2005l) . 



If this interpretation is true, then the shape of the 
bump in the afterglow light curve requires a limb- 
brightened afterglow image, in agreement with the - 
oretical expectations ijGaudi. Granot fc LoeblEoOlj) . 

The size of the afterglow image at a sin- 
gle epoch can be estimated from the quenching 
of diffractive scintillations in the radio afterglow 
llGoodmanll997UFrail et al.ll997tlTavlor et al.ll997t 
IWaxman. Kulkarni fc Frail ll99S|) . The flux be- 
low the self-absorption frequency can also be used 
to constrain the size of the emitting region (e.g., 
Katz fc Piranl Il997t IGranot. Ramirez-Ruiz fc Loebl 
2005). A more direct measurement of the image 
size, as well as its temporal evolution, may be ob- 
tained through very large base-line interferometry in 
the radio (i.e. with the VLB A). This was possible 
for only for one radio afte rglow so far (GRB 030329; 
iTavlor et alj \2004. 120051) . since it requires a rela- 
tively nearby event (z < 0.2) with a bright ra- 
dio afterglow. Nevert heless, it already provides 
interesting constraints (| Oren. N akar fc Piranl l2004t 
IGranot. Ramirez-Ruiz fc Loebll2005L see Figure ITT )l . 
and better observations in the future may help pin 
down the jet structure and dynamics, as well as the 
external density profile. 
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2.6. What causes the Jet Break? 

The jet break in the afterglow light curve has 
been argued to be the combination of (i) the edge of 
the jet becoming visible, and (ii) fast lateral spread- 
ing. Both effects are expected to take place around 
the same time, when the Lorentz factor, F, of the 
jet drops below the inverse of its initial half-opening 
angle, 6q. This can be understood as follows. 

When T drops below 9^ 1 the edge of the jet be- 
comes visible, since relativistic beaming limits the 
region from which a significant fraction of the emit- 
ted radiation reaches the observer to within an angle 
of ~ r _1 around the line of sight (9 < T -1 ). Once 
the edge of the jet becomes visible, then if there is 
no significant lateral spreading, only a small fraction 
(T9j) 2 < 1 of the visible region is occupied by the 
jet, and therefore there would be "missing" contri- 
butions to the observed flux compared to a spheri- 
cal flow. This would cause a steepening in the light 
curve, i.e. a jet break, where the temporal decay in- 
dex asymptotically increases by An = (3— fc)/(4— k). 

When r drops below Oq 1 , the center of the jet 
comes into causal contact with its edge, and the 
jet can in principal start to expand sideways sig- 
nificantly. It has been argued that at this stage it 
would indeed start to expand sideways rapidly, at 
close to the speed of light in its own rest frame. In 
this case, during the rapid lateral expansion phase 
the jet opening angle grows as 9j ~ F -1 and ex- 
ponentially with radius (sec t!2.1|l . This causes the 
energy per solid angle, e, in the jet to drop with 
observed time, and the Lorentz factor to decrease 
faster as a function of the observed time, which re- 
sult a steepening in the afterglow light curve com- 
pared to a spherical flow (where e remains constant 
and r decreases more slowly with the observed time). 
However, in this case a good part of the visible re- 
gion remains occupied by the jet (since TQj remains 
~ 1), so that the first cause for the jet break (the 
edge of the jet becoming visible, and the "missing" 
contributions from outside the edge of the jet) is no 
longer important. Therefore, for fast lateral spread- 
ing, the jet break is caused predominantly since the 
energy per solid angle e decreases with time, and the 
Lorentz factor decreases with observed time faster 
than for a spherical flow. 

It is important to keep in mind, however, that 
numerical studies show that the lateral spreading of 
the jet is very modest as long as it is relativistic (see 
H2.2I and fl2.3|) . This implies that lateral spreading 
cannot play an important role in the jet break, and 
the predominant cause of the jet break is the "miss- 
ing" contribution from outside of the jet, once its 
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Fig. 12. The temporal decay index a as a function of the 
observed time (in days) across the jet break in the light 
curve, for p = 2.5. Upper panel: results in the spectral 
range v > max(^ m , v c ) using a sem i-analytic model with 
no lateral spreading l|Granotl2005r) . for a uniform (k = 0, 
n-cxt = 1 cm -3 ) and wind (k = 2, A, = 1) external den- 
sity profile, with 0q = 0.1 and i?k,iso = 2 x 10 53 erg. 
Lower panel: results for the spectral range v m <v<v a , 
for 8o = 0.2 and a uniform density (k = 0, n cx t = 1 cm -3 , 
-Ek.iso = 10 52 erg); co mpares the result of a semi-analytic 
model |Gnmo1j|2005 ) to th ose of a hydrodynamic simu- 
lation iGranot et al l200lf) . In both panels the dashed 
lines show the asymptotic values of a before and after 
the jet break, for a uniform jet with no lateral spreading, 
for which Aa = (3 - k)/(4 - k). 

edge becomes visible. 

A potential problem with this picture is that if 
the jet half-opening angle remains roughly constant, 
9j « 9q, the asymptotic change in the temporal de- 
cay index is only Aa = 3/4 for a uniform exter- 
nal medium (k = 0) or even smaller for a wind 
(Aa = 1/2 for k = 2), while the values inferred 
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from obs ervations are in most case s larger (see Fig- 
ure 3 of IZeh. Klose fc KannlhoOfift . This apparent 
discrepancy may be reconciled as follows. While the 
asymptotic steepening is indeed Aa = (3 — fc)/(4— fc) 
when lateral expansion is negligible, the value of the 
temporal decay index a (where F v oc t~ a ) initially 
overshoots its asymptotic value. Since the temporal 
baseline that is used in order to measure the post-jet 
break temporal decay index ct2 is typically no more 
than a factor of several in time after the jet break 
time 8 , tjet, the value of a during this time is larger 
than its asymptotic value «2- This causes the value 
of Aa that is inferred from observations to be larger 
than its asymptotic value. 

The overshoot in the value of a just after the 
jet break time can nicely be seen in Figure ^1 and 
is much more pronounced in the light curves calcu- 
lated using the jet dynamics from a hydrodynamic 
simulation, compared to the results of a simple semi- 
analytic model. The cause of this overshoot is that 
the afterglow image is limb-brightened (see Figure 
EjJ) and therefore the outer edges of the image which 
are the brightest are the first region whose contri- 
bution to the observed flux is "missed" as the edge 
of the jet becomes visible. The overshoot is larger 
the more limb-brightened the afterglow image (e.g., 
for v > max(i/ m ,i/ c ) in the upper panel of Figure 
1121 compared to v m < v < v c in the lower panel of 
Figure El- For a wind density (fc = 2) the limb- 
brightening is smaller compared to a uniform den- 
sity (fc = 0), at the same power law segment of the 
spectrum (see Figure EJ), and the Lorentz factor T 
decreases more slowly with the observed time. Be- 
cause of this no overshoot is seen in the semi-analytic 
model shown in the upper panel of Figure 1121 for a 
wind density profile (k = 2), and the jet break is 
smoother and extends over a larger factor in time. 
The asymptotic post-jet break value of the temporal 
decay index (a^) is approached only when the visi- 
ble part of the afterglow image covers the relatively 
uniform central part, and not its brighter outer edge. 

The jet break in light curves calculated from 
hydrodynamic simulations is sharper than in semi- 
analytic models (where the emission is taken to be 
from a 2D surface - usually a section of a sphere 
within a cone). In semi-analytic models the jet break 
is sharpest with no lateral expansion, and becomes 
more gradual the faster the assumed lateral expan- 
sion. For example, in the lower panel of Figure 1121 
where the viewing angle is along the jet axis and the 



external density is uniform, most of the change in 
the temporal decay index a occurs over a factor of 
~ 2 in time for the numerical simulation, and over 
a factor of ~ 3 in time for the semi-analytic model 
(which assumes no lateral expansion; the jet break 
would be more gradual with lateral expansion). For 
both types of models, the jet break is more grad- 
ual and occurs at a somewhat later time for viewing 
angles further away from the jet symmetry axis but 
still within its initial opening angle, although this ef- 
fect is s omewhat more pronounced in semi-a nalytic 
models l|Granot et al.ll200ll iRossi et al.ll2004|) . 



3. THE JET STRUCTURE 

Since the initial discovery of GR B after- 

J lows in the X-ray llCosta et al.l Il997ft. optical 
van Paradiis et all Il997ft . and radio l|Frail et alJ 
Il997ft . many afterglows have been detected and the 
quality of individual after glow light curves h as im- 
proved dramatically (e.g.. iLipkin et ai1l2004|) . De- 
spite all the observational and theoretical progress, 
the structure of GRB jets remains largely an open 
question. This question is of great importance and 
interest, since it is related to issues that are funda- 
mental for our understanding of GRBs, such as their 
event rate, total energy, and the requirements from 
the compact source that accelerates and collimatcs 
these jets. 

In H3.ll a brief overview is given of the main jet 
structures that have been discussed in the literature 
and the motivation for them. This is followed by a 
discussion of the different methods that have been 
applied for constraining the jet structure from ob- 
servations, which include statistical studies f H3.2|) . as 
well as the evolution of the linear polarization of the 
afterglow emission f t!3.3|) . and the shape of the after- 
glow light curves ( ^3.4[l . The afterglow light curves 
from viewing angles outside the initial jet aperture 
are discussed in H3. 51 along with possible implications 
for X-ray flashes and for the jet structure, while H3.6I 
briefly mentions the search for orphan afterglows. 
Some implications of recent Swift observations are 
discussed in iJ3.7l 

3.1. Existing Models for the Jet Structure 

The leading models for the jet s tructure 
are (i) the unif orm jet (UJ) model jRhoadd 
19971 Il999t IPanaitescu fc Meszarod Il999t 



8 This is usually because the flux becomes too dim to de- 
tect above the host galaxy, or since a supernova component 
becomes dominant in the optical, etc. 



Sari. Piran fc HalperrJ Il999t I Kumar fc Panaitescul 



200C : iModerski. Sikora fc Bulikll2000t iGranot et al 
2001112002ft . where the energy per solid angle, e, and 
the initial Lorentz factor, Tq, are uniform within 
some finite half-opening angle, 6j, and sharply drop 
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outside of 9j ; and (ii) the universal structured jet 
(USJ) model llLipunov. Postnov fc ProkhorovlboOlt 
iRossi et all l2002i IZhang fc Meszarosl l2002h . where 
e and To vary smoothly with the angle 9 from the 
jet symmetry axis. In the UJ model the different 
values of the jet break time, tj , in the afterglow 
light curve arise mainly due to different Oj (and to 
a lesser extent due to different ambient densities). 
In the USJ model, all GRB jets are intrinsically 
identical, and the different values of tj arise mainly 
due to different viewing angles, 6* b s , from the jet 
axis. 



Th e observed correlation, tj o c E } (|Frail et all 

l200lt iBloom. Frail fc Kulkarnl 120031) . implies a 
roughly constant true energy, E, between differ- 
ent GRB jets in the UJ model, and e oc 9~~ 2 out- 
side of some cor e angle, r , in the USJ model 
(|Rossi et al.ll2002t IZhang fc Meszarosl l2002j) . This 
is assuming a constant efficiency, e 7 , for producing 
the observed prompt gamma-ray (or X-ray) emis- 
sion. If the efficiency depends on 9 in the USJ 
model, for example, then different power laws of 
e wit h 9 are possible (|Guetta. Granot fc Begelmanl 
l2005lh such as a core with wings where e oc 



as is 



o btained in simulations of the col- 
lapsar model ilZhang. Woosley fc MacFadvenl l 2003t 



i/Vooslev 

IZhang. Wooslev fc Hegerll2004 . 

The jet structure was initially envisioned to be 
uniform since this is the simplest jet structure, and 
arguably also the most natural. Furthermore, it 
also p redicted a jet break in the afterglow light 
curve jRhoadsl ITM ISari. Piran fc tialneml 

Il999|), which was soon thereafter confirmed observa- 
tionallv dFruchter et al.lll999t iKulkarni et"aDll999ai 
IStanek et al.lll999|) . The original motivation for the 
USJ model was the conceptual simplicity of a uni- 
versal intrinsic structure for all GRB jets, where the 
observed differences (namely in the jet break times) 
arise due to different viewing angles (instead of being 
attributed to an intrinsic difference - in the jet half- 
opening angle - as in the UJ model) . Its exact struc- 
ture was motivated by the requirement to reproduce 
the observed afterglow light curves and correlations 
with the prompt GRB emission. It had later been 
suggested on theoretical grounds that a jet struc- 
ture with a narrow core and wings where e oc 9~ 2 is 
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9 In fact, the expression for tj is similar to that for a uniform 
jet with e — > e(6 = (9 b s ) and 6j — i 9 ^, s 

10 Simulations of jets launched by an accretion torus - black 
hole system, where the jet does not propagate through a 
stellar envelope, as is expected to arise in binary merger 
scenarios which might be relevant to GRBs of the short- 
hard class, produce a rough ly uniform jet with sharp edges 
lAlov. Janka fc Kullerll2005l) . 
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Fig. 13. An illustration of various jet structures that have 
been discussed in the literature, in terms of the distribu- 
tion of their energy per solid angle, e = dE/dQ., with 
the angle 6 from the jet symmetry axis, both in semi- 
logarithmic scale [main figure) and in log- log scale {big 
inset). Both the normalization of dE/dQ and the typical 
angular scale may vary in most models, and their val- 
ues s hown here wer e chosen to be more or less "typical" . 
(from lGranotlbOOED . 

expected i n highly magnetized Poynting flux domi- 
nated jets l|Lvutikov fc Blandfordl2002ll2003|) as well 
as in low magnetization hydr o-dynamic jets in the 
conte xt of the collaos ar model i La^za^fc^JBe^ehnan 
2005 : see, however, iMorsonv. Lazzati fc Begelman 



2006) 



Other jet structures have also been proposed in 
the literature. Figure ^] illustrates different jet 
structures that have been discussed in the litera- 
ture in terms of their distribut ion of e(9). A jet 
with a Gaussian angular prof ile ijZhang fc Meszarosl 
120021: iKumar fc Granotll2003|) may be thought of as 
a more realistic version of a uniform jet, where the 
edges arc smooth rather than sharp. A Gaussian 
e(9) oc exp(— 9 2 /29 2 ) is approximately intermediate 
between the UJ and USJ models, but it is closer to 
the UJ model than to the USJ model with e oc 9~ 2 
in the sense that for a Gaussian e(9) the energy in 
the wings of the jet is much smaller than in its core, 
whereas for a USJ with e oc 9~ 2 wings there is equal 
energy per decade in the wings, and therefore the 
wings contain more energy than the core (by about 
an order of magnitude). 

Another jet structure that received some at- 
tention recently is a two-component jet model 
JPedersen et alJll998t lFrail et al-lhoOOl iBerger et al. 
2003bl: iHuang et all [2004t IPeng. K5nigl fc Gra~ 
20051 IWu et all 12005 ) with a narrow uniform jet 



of initial Lorentz factor Tq > 100 surrounded by a 
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wider uniform jet with r ~ 10 — 30. Theoretical 
motivation for such a jet structure has been found 
both in the context of the cocoon in the c ollapsar 
model ( Ramirez-Ruiz. Celotti fc R.eesl l2002j) and in 
the con text of a hydromagnetically dri ven neutron- 
rich jet l|Vlahakis. Peng fe K5nig|2003|) . This model 
has been invoked in order to account for sharp 
bumps (i.e., fast rebrightening episod es) in the af- 
terglow light curves of G RB 030329 feerger et all 
l2003blt and XRF 030723 llHuang et all 1200^ . A 
different motivation for proposing this jet structure 
is in order to account for the energetics of GRBs 
and X-ray flashes and reduce the high efficiency re- 
quirements from the prompt gamma-ray emission 
(IPeng. Konigl fc GranotlboOfil) . 

More "exotic" jet structures have also been con- 
sidered. One example is a jet with a cross sec- 
tion in the shape of a "ring," sometimes referred t o 
as a "hollow cone" jLevinson fc Eichlerlll993[ EoO 
Eichler fc Levinson1l2003l l2004t lLazzati fc Begelma 



20051) . which is uniform within 9 C < 9 < 9 C + AO 
where A9 <C 9 C . An other example is a "fan"- 
or "sheet" -shaped jet ((Thompson! I2005T) where a 
magnetoccntrifugally launched wind from the proto- 
neutron star, formed during the supernova explosion 
in the massive star progenitor, becomes relativistic 
as the density in its immediate vicinity drops and is 
envisioned to form a thin sheath of relativistic out- 
flow that is somehow able to penetrate through the 
progenitor star along the rotational equator, forming 
a relativistic outflow within A9 <C 1 around 9 = n/2 
(or 9 C = 7r/2- A9/2). 11 

3.2. Statistical Studies 

One approach for constraining the jet structure 
is through statistical studies of the prompt emis- 
sion. In particular, a convenient observable is the 
log N — log 5* distribution, where N is the number of 
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to assume 



both the intrinsic GRB event rate (which is usually 
assumed to follow the star formation rate), and the 
luminosity function which depends on the jet struc- 
ture through the isotropic equivalent luminosity L(9) 
as a function of the angle 9 from the jet symmetry 
axis. The latter depends in turn on the angular dis- 
tribution of the energy per solid angle in the jet e(9) 
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11 This has been suggested as a possible jet structure within 
this model, but the final jet structure is by no means clear, 
and other jet structures might also be possible within this 
model (T. A. Thompson 2005, private communication). 



0.1 0.2 0.3 0.4 0.5 

6 (rad) 



Fig. 14. The probability distribution, dN/dzdln8, of 
the observed GRB rate as a function of redshift z and 
viewing angle 8, as predicted by the USJ model (from 
iNakar. Granot fc Guettall2004l) . The white contour lines 
confine the minimal area that contains 1 a of the total 
probability. The circles denote the 16 GRBs with known 
z and 6 from the sample of iBloom. Frail fc Kulkarnll 
ll2003l). (a) The mode l para meters are similar to those 
of iPerna Sari fc Fraill (120031) . This figure is the 2D re- 
alization of their Figure 1. (b) Here a limited range in 
redshift is used, 0.8 < z < 1.7 (containing 10 out of the 
16 data points), in order to minimize redshift selection 
effects and reduce the sensitivity of the results to the 
unknown GRB rate. Measurement errors of 20% in ln# 
(o'lnB = 0.2) were included and a log-normal distribution 
in the effective duration that they deduced from obser- 
vations. 



and the gamma-ray efficiency e 7 (0), whose product 
(times 47r) provides the isotropic equivalent energy 
output in gamma-rays, Ey t - mo (6), and the assumed 
distribution of the peak isotropic equivalent lumi- 
nosity for a given Ey[ S0 . 

For the USJ model e(9) and e 7 (0) are assumed 
to be power laws in 9 and one needs to specify their 
power law indexes, as well as the jet core angle 9 r 
and outer edge # max l|Guetta. Granot fc Begelmanl 
120051) . For the UJ model one needs to specify 
the distribution of jet half-opening angles 9 , which 
can be taken to be a power law distribution in 
the range 9 mm < 6q < # max . These simple forms 
of the USJ and UJ models are degenerate, since 
they both produce a power law luminosity function 
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l|Guetta. Granot fc Begelmanll2005() . which provides 
an adequate fit to the data. The total rate of en- 
ergy release in gamma-rays must be the same in 
both models (and match the observed rate), where 
in the USJ model it is released in fewer more ener- 
getic events (by a factor of ~ 10; s ee equation 14 of 
iGuetta. Granot fc Begelmanll2005|) . while in the UJ 
model it is released in more numerous less energetic 
events (i.e. the UJ model predicts a larger intrinsic 
event rate). 

Another statistical approach for constraining the 
jet structure is through the distribution of the ob- 
served number of GRBs A as a function of the angle 
9 that is inferred from the observed jet break times in 
the afterglow light curves, where 9 corresponds to the 
jet half-opening angle 9$ in the UJ model and to the 
viewing angle # bs hi the USJ model. The observed 
dN/d9 distribution agrees reaso nably well with the 
predi ctions of the USJ model l|Perna. Sari fc Fraill 
120031) . which had been argued to support the USJ 
model (since in the competing UJ model there is 
an additional freedom in the choice of the proba- 
bility distribution for 9q which would make it eas- 
ier to fit these observations). However, when the 
known redshifts z of the GRBs in the same sam- 
ple are also taken into account, then the predictions 
of the USJ model for the two dimensional distribu- 
tion of observed GRBs with 9 and z, dN/dzd9, is 
found to be in very poor agreement with observa- 
tions (|Nakar. Granot fc Guettall2004|) . This can be 
best seen for a relatively narrow range in z (see lower 
panel of Figurc H4*|) , in which the USJ model predicts 
that most GRBs should be near the upper end of 
the observed range in 9, while in the observed sam- 
ple most GRBs are near the lower end of that range. 
Since the available sample was very inhomogeneous 
(i.e., involved many different detectors), it should be 
taken with care and cannot be used to rule out the 
USJ model. Nevertheless, it strongly disfavors the 
USJ model. 

3.3. Linear Polarization 
Linear polarization at the level of a few percent 
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IWiiers et al.lll999t 

as is illustrated in 
Figure 1151 This was considered as a confirmation 
that synchrotron radiation is the dominant emission 
mechanism in the afterglow. The most popular ex- 
planation for the observed linear polari zation had 
been synchrotron emissio n from a jet l)Saril Il999t 
iGhisellini fc Lazzatilll999l) . In this model the mag- 
netic field is produced at the afterglow shock and 



possesses axial symmetry about the shock normal. 
In this picture there would be no net polarization 
for a spherical outflow, since the polarization from 
the different parts of the afterglow image would can- 
cel out, and a jet geometry together with a line of 
sight that is not along the jet axis (but still within 
the jet aperture, in order to see the prompt GRB) is 
needed in order to break the symmetry of the after- 
glow image around our line of sight. 

For a uniform jet (the UJ model) this predicts 
two peaks in the polarization light curve around the 
jet bre ak time tj, if T9j < 1 decreases with time a t 



t > tj l|Ghisellini fc Lazzatilll999llRossi etal 



2004) 



1999) 



or even three peaks if T9j ~ 1 at t > tj jSarl 
where in both cases the polarization vanishes and 
reappears rotated by 90° between adjacent peaks. 
The latter is a distinct signature of this model. For a 
structured jet (the USJ model), the polarization po- 
sition angle is expected to remain constant in time, 
while the degr ee of polarization peaks near the jet 
break time tj IjRossi et al.ll2004l) . A similar qualita- 
tive behavior is also expected for a Gaussian jet, or 
other jet structures with a bright core and dimmer 
wings (although there are obviously some quantita- 
tive differences). 

The different predictions for the afterglow polar- 
ization light curves for different jet structures raise 
the hopes that afterglow polarization observations 
may constrain the jet structure. In practice, how- 
ever, the situation is much more complicated, mainly 
since the observed polarization depends not only on 
the jet geometry, but also on the magnetic field con- 
figuration in the emitting region, which is not known 
very well. For example, an ordered magnetic field 
component in the emitting region (e.g. due to a small 
ordered magnetic field in the external medium) may 
dominate the polarized flux, and therefore the po- 
larization light curves, even if it is sub-dominant in 
the emitting region compared to a random (shock 
generated) magnetic field compon ent in terms of the 
total energy in the magnetic field I Granot fc KonigJ 
120031) . Other models for afterglow polarization in- 
clude a magnetic field that is coherent over patches 
of a si ze comparable to that of cau sally connected re- 
gions (|Gruzinov fc Waxmanlll999l). and polarization 
that is induced by microlensing JLoeb fc Pernall998l) 
or by scintillations in the radio ))]VieTlve^h^fcLoebl 
119991) . 

An additional complication arises in afterglow 
light curves that exhibit variability, since a vari- 
able afterglow light curve is expected to be ac- 
companied by a variable polarization light curve, 
both in the degree of polarization and in its 
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Fig. 15. Summary of afterglow p olarization measurement 
up to 2002 rfrom lCovinoll2004l) . Top panel: degree of 
polarization and position angle for all the positive detec- 
tions, i.e. upper limits are excluded. Bottom panel: Q 
and U Stokes parameters for all the available data, i.e. 
including upper limits. 



position angle l|Granot fc Konid |2003|) . This 
prediction has been confirmed in GRB 021004 
(|Rol et alJ l2003[) . where it had been interpreted 
both in the context of angular i nhomogeneities 
withi n the jet (i.e. a "patchy shell" , iNakar fc Orenl 
|2004|) and as discrete episodes of energy injec- 
tion into the afterglow shock (i.e. "refreshed 
shocks" , iBiornsson. Gudmundsson fc Johannessor 



2004J), and later also in GRB 030329 (|Greiner et al 
200 



Perhaps the best monitored polarization light 
curve of a smooth afterglow, which does not suf- 
fer from the complications mentio ned above, is 
GRB 020813 (|Gorosabel et al.ll2004|) where the po- 
larization position angle is roughly constant while 
the degree of polarization decreased by a factor of 
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Fig. 16. Fits of the predicted polarization for different 
jet structures and magnetic field configuration to the op- 
tical afterglow pola rization data for GRB 020813 (from 
lLazzati et alJl2f)04) . The upper panel shows the degree 
of polarization and position angle around the jet break 
time, which is also where the observations are concen- 
trated, while the bottom panel shows the degree of polar- 
ization over a wider range in time and demonstrates that 
the differences between the various models are more pro- 
nounced at t < tj . The shaded region denotes the allowed 
range for the jet break time tj. The external density is 
taken to be either uniform (ISM, k = 0) or a stellar 
wind (Wind, k = 2). The different models are labeled 
by a three letter acronym where the first letter describes 
the magnetic field ('H' is for hydrodynamic, i.e. sock 
produced magnetic field, while 'M' is for magnetized, i.e. 
with an ordered toroidal magnetic field) while the second 
letter describes the jet structure ('S' for a structured jet, 
and 'H' for a homogeneous or uniform jet). 



JETS IN GAMMA-RAY BURSTS 



17 



~ 2 from ~ 0.5tj to ~ 2tj. The constant position 
angle across the jet break time tj disfavors a uni- 
form jet with a shock generated magnetic field l|Saril 
Il99flt iGhisellini fc Lazzatilll999l), as well as patches 
of uniform field l|Gruzinov fc Waxmax]ll999|) where 
the position angle (as well as the degree of polariza- 
tion) is expected to change randomly on time scales 
At <t 

lLazzati et al.l (|2004|) have contrasted different 
models for the jet structure and magnetic field con- 
figuration with the polarization data for GRB 020813 
(see Figure lTT)|) . and concluded that the data support 
either (i) a structured jet (USJ) or a jet structure 
where most of the jet energy is in a narrow core 
while its wings contain less energy (such as a Gaus- 
sian jet) with a shock produced magnetic field (where 
the field is not purely in the plane of the shock but 
still possesses significant anisotropy) or (ii) a uni- 
form jet 12 with an ordered toroidal magnetic field 
component which dominates the polarized flux to- 
gether with a random magnetic field component that 
dominates the total flux (and the total magnetic en- 
ergy) in order for the polarization not to exceed the 
observed value. 

We conclude that while the afterglow polariza- 
tion light curves may provide useful constraints on 
the jet structure and the magnetic field configuration 
in the emitting region, it is in practice rather difficult 
to constrain each of these two ingredients separately. 
That is, in order to obtain tight constraints on the jet 
structure, strong assumptions must be made on the 
magnetic field configuration, and vice versa. Nev- 
ertheless, as discussed above, interesting constraints 
have already been derived from existing data. 

3.4. Afterglow Light Curves 

The shape of the afterglow light curves is an 
important and relatively robust diagnostic tool for 
constraining the jet structure. The afterglow light 
curves (at least starting from a few hours after the 
GRB) are typically described by an initial power law 
flux decay F v oc t~ a with 0.7 < ot\ < 1.5 which steep- 
ens into a sharper power law decay (1.6 < a-> > < 2.8 ) 
at the jet break time tj ijZeh. Klose fc KannlfeOOtj) . 
Figure 1171 shows an example of the very well mon- 
itored jet break of GRB 030329. The jet break in 
the light curve is usually rather sharp (most of the 
steepening occurs within a factor of a few in time) 

^Lazzati et al.l J2004T) find that a structured jet (USJ) with 
an ordered toroidal field component that dominates the po- 
larization is still hard to rule out, even though it provides a 
poor fit to the polarization data of GRB 020813, due to the 
model uncertainties regarding the mixing of such an ordered 
field component with the shock generated random field. 
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Fig. 17. The jet break of GRB 030329 at three dif- 
ferent optical bands, with e xquisite temporal sampling 
(from iGorosabel et al.l l200et) . The upper panel shows 
a wider time span complemented by data from other 
groups, while t h e mi ddle panel shows the data from 
IGorosabel et al.l (1200(1) . and the lower panel shows the 
residuals relative to a fit to a smooth jet break model. 
The steepening in the temporal decay during the jet 
break is smooth and achromatic. 

and the increase in the temporal decay index, Aa = 
ct2 — eg , is typically in th e range 0.7 < Aa < 1.4 
l|Zeh. Klose fc Kannll2006|) . Different jet structures 
may be tested by their ability to reproduce these ob- 
served properties. Below we describe the resulting 
constraints for various jet models. 

Uniform Jet: Figure ^| shows the afterglow 
light curves for an initially uniform jet whose evo- 
lutio n is calculated usin g a hydrodynamic simula- 
tion l|Granot et al-lEoOlf) . The initial conditions arc 
a cone of half-opening angle 8q taken out of the 
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Fig. 18. Afterglow light curves for an initially uniform jet 
whose evolution is calculated using a hydrodynamic sim- 
ulation, at different observe frequ encies, for an observ er 
along the jet symmetry axis (from iGranot et aT1l200ll) . 



spherical self-similar iBlandford fc McKeel ( 197(t) so- 
lution (see i |2.3|l . The shape of the light curves is 
nicely consistent with those observed in real after- 
glows, particularly in terms of the sharpness of the 
jet break, where the observed diversity can be at- 
tributed to different viewing angles within the jet 
aperture, 6> b s < #o ( see upper panel of Figure 

Gaussian Jet: the afterglow light curves for a 
Gaussian jet where e(9) = en exp(— 9 2 /29%), that is 
observed at viewing angles inside the core of the jet 
($obs < Or.) are rather similar to those for a uni- 
form iet llKumar fc Gra"notll2003l iRossi et allEool 
IGranot. Ramirez-Ruiz fc Pernal 120051) . and reason- 
ably consistent with afterglow observations. 

Structured Jet: one can consider a jet with a 
narrow core and wings where both e and T are ini- 
tially power laws in the angle 9 from the jet symme- 
try axis: e oc 9~ a and To oc 9~ b , outside some (nar- 
row) core angle 9 C . A comparison between the result- 
ing afterglow light curves and afterglow observations 
can then be used to constra in the power law indexes 
a and 6. Such an analysis <)Granot fc Kumarll2003h 
suggests that < b < 1 and a « 2 (or 1.5 < a < 2.5, 
see Figure 1*031 . 

The upper limit on b comes from the fact that a 
large value of b implies a small initial Lorentz fac- 
tor at large viewing angles, Tr>(9 h s 3> &c), since it 
is hard for Tq(9 = 0) to exceed 10 4 . A small initial 
Lorentz factor along the line of sight at large view- 
ing angles would result in a large decelerations time 
along the line of sight and therefore an initially rising 
light curve, up to relatively late times, which is not 
seen in observations. Furthermore rn(6> b s ) > 100 is 
needed in order to produce the prompt gamma-ray 
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Fig. 19. Afterglow light curves for a structured jet 
where initially e oc 8~ a and T oc 8~ b ou tside some 
core ane le 8 C (from IGranot fc Kumar! 120031) . The dif- 
ferent curves, from top to bottom, are for viewing an- 
gles 8 obs = 0.01, 0.03, 0.05, 0.1, 0.2, 0.3, 0.5, where 8 C = 
0.02, p = 2.5, e e = e B = 0.1, n cxt = 1 cm" 3 , T(8 = 
0, to) = 10 3 , and the total energy in the jet is 10 52 erg. 
In model 1 e(8) does not change with time. In model 2, 
e(8, t) evolves such that it is proportional to the average 
over its initial distribution, e(8,to), over the range in 8 
out to which a sound wave coul d propagate from to up 
to t (see IGranot fc Kumarll2003l . for details). 

emission. Figure ED shows light curves for 6 = 
and a = 1, 2, 3 and different viewing angles. For 
a = 1 the change in the temporal decay index across 
the jet break, Act, is too small (compared to its ob- 
served values) and the post-jet break decay slope is 
not steep enough. For a = 3 there is either a very 
pronounced flattening in the light curve before the 
jet break time or the temporal decay slope after the 
jet break is extremely steep (neither of which is seen 
in afterglow observations). This suggests a w 2 (or 
1.5 < a < 2.5). 

As can be seen from Figure [201 even for a = 2 
and 6 = there is a flattening in the light curve 
before the jet break time, which becomes more pro- 
nounced at large viewing angles (0 o bs 3> C ). This 
arises since the bright core of the jet becomes visible, 
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Fig. 20. Optical (R-band) light curves (upper panel) and 
the temporal decay index (lower panel) for a struct ured 
jet with no lateral spreading (from lRossi et af]|2004f) . for 
a = 2, b = 0, £k,iso(0 = 0) = 2 x 10 54 erg, n ext = 1 cm" 3 , 
e e = 0.01, e B = 0.005, 0c = 1°- The different curves are 
for different viewing angles, # bs, which are labeled in 
units of the jet core angle, 9 C . 



while the value of e along the line of sight is much 
smaller in comparison for large viewing angles, which 
more than compensates for the relatively small frac- 
tion of the visible region that is occupied by the jet 
core. That is, the average energy per solid angle 
in the observed region initially increases with time 
until the core of the jet becomes visible, and then 
decreases with time. Such a flattening has not been 
observed so for, but since it is most pronounced at 
large viewing angles for which the jet break time is 
large and the flux around that time is low, the ob- 
servations in this parameter range are usually rather 
sparse, making it hard to rule out this model on these 
grounds. 

Two Component Jet Model: the light curves 
for this j et structure have been calcu lated an- 
alytically llPeng. Konigl fc Granoti 120051) or scmi- 
analytically llHuang et all2004IWu et al.l2005h . and 

it has been suggested that this model can account for 
sharp bumps (i.e., fast rcbrightcning ep isodes) in the 
aftergl ow light curves of GRB 030329 feerger et al.1 
l2003bh and XRF 030723 llHuang et al.ll2004l) . It has 
been later demonstrated, however, that effects such 
as the modest degree of lateral expansion that is ex- 
pected in impulsive relativistic jets and the gradual 
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Fig. 21. Light curves for a jet with an angular struc- 
ture of a ring, where e is uniform for 8 C < 6 < C + A9, 
and sharply drops outside of this range in 9, for various 
fractional widths, viewed from within the jet. The upper 
line is for a uniform jet viewed from along its symme- 
try axis (9 C = ^obs = 0, AO — 0.2) and is included for 
comparison, while the other lines are for a ring shaped 
jet with 9 C = 0.1 and 6 C /A9 = 1, 2, 3, 5, 10, from top to 
bottom, viewed from 9 h s = C + A9/2. A uniform ex- 
ternal density (k = 0), p = 2.5, a nd no lateral expansion 
are assumed (from lGranotll2005ll . 

hydrodynamic transition at the deceleration epoch 
smoothen the resulting features in the afterglow light 
curve, so that they cannot produce featu res as sharp 
as those mentioned above l|Granotll2005|) . One of the 
motivations for a two component jet was to alleviate 
the efficiency requirements on the prompt g amma- 
ray emission llPeng. Konigl fc Granotl boOSI) if the 
wide jet dominates the total energy and late time af- 
terglow emission, while the narrow jet is responsible 
for the prompt gamma-ray emission. However, the 
more recent Swift observations, which show a flat de- 
cay p hase in the early X-ray afterglow l|Nousek et alJ 
l20f)d) . are inconsistent which this picture and sug- 
gest a high gamma-ray efficiency also for this jet 
model under the standard assumption s of afterglow 
theory l|Granot. K5nigl fc Piranll2006^ . 

Ring Shaped Jet: Figurcl2*Tlshows light curves 
for a ring shaped jet, viewed from within the cmit- 
ting region, fo r different fractional width of the ring 
llGranod E)?)! . For a thin ring, whose width is 
smaller than its radius, A9 < 8 C , the jet break is 
divided into two distinct and smaller breaks, the 
first occurring when TA9 ~ 1 — 2 (i.e. when the 
width of the jet becomes visible) and the second 
when T6 C ~ 1/2 (i.e. when the whole jet becomes 
visible). This is inconsistent with the large steepen- 
ing that is observed across the jet break, and sug- 
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Fig. 22. Light curves for a jet with an angular structure 
of a thin fan, with an opening angle of Af9 = 0.1 centered 
on 6 = 71-/2, i-e. C = tt/2 - A0/2 = tt/2 - 0.05 (from 
lOranodEfflsh . 



gests a "thick" ring, with AO > C . Even for such 
a thick ring with AO = C the jet break is still not 
quite sharp enough to match observations according 
to the model shown in Figurc l2*Tl however, the inclu- 
sion of some lateral expansion which would tend to 
smoothen the edges of the jet and fill in the center 
of the thick ring would probably be enough to make 
the light curve for A6> > C consistent with observed 
afterglow light curves. The fact that the beaming 
cone extends out to an angle of ~ from the edge 
of the jet further contributes to smoothing the jet 
break for a thick ring jet, and contributes to the ob- 
ser ved flux from lines of sight near the edge of the 
jet (jEichlerll2005^ . 

Fan Shaped Jet: Figure [221 shows the light 
curves for a jet in the shape of a thin fan, for different 
viewing angles. This jet structure is a limiting case 
of the ring shaped jet where C = tt/2 — AO/2, and 
the jet occupies an angle of AO <C 1 centered around 
= 7r/2. In this case the second jet break from the 
thin ring jet that had been discussed above occurs 
only when the jet becomes sub-relativistic, and is in- 
distinguishable from the non-relativistic transition. 
This leaves only one jet break, when rA6> ^1 — 2 
and the width of the jet becomes visible. The steep- 
ening in the flux decay rate across this jet break is at 
most half of that for a uniform jet (exactly half with 
no lateral expansion, and slightly less than half with 
rapid lateral ex pansion), and is too small compared 
to observations l)Granoti2005|) . This practically rules 
out a jet in the shape of a thin fan. 



3.5. Off-Axis Viewing Angles 
When the jet has relatively sharp edges, 
then the observed emission from viewing an- 
gles outside of the jet aperture is very dif- 
ferent than that f rom viewing angles within 
the jet ap erture llPanaitescu. Meszaros fc Reesl 
19981 IPanaitescu fc MeszarosT Il999t 

Moderski. Sikora fc Bulild l200fi iGranot et all 
200ll l2002t iDalal. Griest fc Pruitl |2002|) . If the 
line of sight is at an angle of A# obs outside the 
edge of the jet (i.e. from the nearest point in 
the jet with significant emission; for a uniform 
jet A$ bs = #obs — 0j) then at early times when 
rA# bs 3> 1 the emitted radiation is strongly 
beamed away from the observer. This is since 
the emission is roughly isotropic in the comoving 
frame of the emitting material in the jet, and 
it is collimated within an angle of around 
its direction of motion in the lab frame, due to 
relativistic beaming (i.e. aberration of light). As 
the jet decelerates, the beaming of radiation away 
from the line of sight becomes weaker, resulting in 
a rising flux. Eventually, when rAf5 b s decreases to 
~ 1, the beaming cone of the jet encompasses the 
line of sight (see Figure I23f) and the observed flux 
peaks and starts to decay, asymptotically joining 
the decaying light curve for viewing angles within 
the initial jet aperture. 

The afterglow light curves for different jet struc- 
tures, dynamics, and viewing angles are shown in 
Figure 1241 For an initially uniform jet with sharp 
edges whose evolution is calculated using a hydro- 
dynamic simulation the rise to the peak in the light 
curve for viewing angles outside the initial jet aper- 
ture (# bs > Oq) is much more gradual compared to 
a semi-analytic model where the jet remains uniform 
with sharp edges and no lateral expansion. 13 This is 
because of the mildly relativistic material at the sides 
of the jet whose emission is not strongly beamed 
away from such off-axis lines of sight (0 o \, s > 0q) at 
early times, and therefore dominates the observed 
flux early on. 

For a Gaussian jet, if both e(0) and Fq(0) have 
a Gaussian profile (corresponding to a constant rest 
mass per solid angle in the outflow), then the af- 
terglow light c urves are rather similar to those for 
a uniform jet ijKumar fc Granotl l2003|) . If, on the 
other hand, e(0) is Gaussian while 14 ^o(O) = const, 
then the light curves for off-axis viewing angles 

13 If lateral expansion is included in such a semi-analytic 
model the rise to the peak in the light curve for (2 b s > 8o 
becomes even steeper, since the beaming cone of the jet ap- 
proaches and eventually engulfs the line of sight faster. 

14 This corresponds to a Gaussian angular distribution of 
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Fig. 23. An illustrative diagram of the emission from a 
uniform relativistic jet with sharp edges and half-opening 
angle 8o, that is seen by an off-axis observer whose line of 
sight makes an angle Sobs > #o with the jet axis. Because 
of relativistic beaming (i.e. aberration of light) the emis- 
sion from each part of the jet is beamed into a narrow 
cone of half-opening angle 7 -1 around its direction of mo- 
tion in the observer frame. During the prompt emission 
(and the very early afterglow) the Lorentz factor of the 
jet is large (7 > 50) and therefore most of the radiation 
is strongly beamed away from the line of sight. In this 
case, the little radiation that is observed comes mainly 
from near the edge of the jet, at the point closest to the 
line of sight. As the jet decelerates 7 decreases with time 
and the beaming cone grows progressively wider, causing 
the radiation to be less strongly beamed, resulting in a 
rising light curve. The light curve peaks when 7 drops 
to ~ (Sobs — do) -1 as the line of sight enters the beaming 
cone of the emitting material at the edge of the jet (the 
middle beaming cone in the figure), and subsequently 
decays with time, asymptotically approaching the light 
curve for an on-axis observer (0 o hs < dp) at later times, 
(from iGranot. Ramirez-Ruiz fc Pernall2005t) . 

(i.e., outside the core of the jet) have a much 
higher flux at early times, compared to a Gaus- 
sian Tq(9) or a uniform jet (see the bottom two 
panels of Figure I24|) , due to a dominant contribu- 
tion from the emitting material along the line of 
sight which has an early deceleration tim e in this 
case ijGranot. Ramirez- Ruiz fe Pernall2005j) . Such a 
jet str ucture was consider ed as a quasi-universal jet 
model l|Zhang et alJl2004(l . 

It has been suggested that a uniform jet with 
sharp edges viewed from slightly outside of its 
edge (60 < 8 bs ^ 20 o ) would result in an X- 
ray flash (XRF) or X-ray rich GRB (XRG RB) 



( Yamazaki. Ioka fe Nakamural2002Ll2003Ll2004() . be- 
cause of the smaller blue-shift of the prompt emission 
compared to viewing angles within the jet (0 o bs < 
9q). This has also been found to nicely explain the 



the rest mass per solid angle, i.e., very little mass near the 
outer edge of the jet, which is the opposite of what might be 
expected due to mixing near the walls of the funnel in the 
massive star progenitor. 
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Hydrodynamic simulation of an 

initially uniform sharp edged jet 




Fig. 24. Afterglow light curves for different jet structures, 
dynamics, and viewing angles. The top panel is from an 
initially uniform jet with sharp edges whose evolution is 
calculated u sing a hydrodynami c simulation (taken from 
Figure 2 of IGranot et al.l l2002l) . The remainin g three 
panels are taken from Figure 5 of IGranot! ll2005tl . where 
a simplified jet dynamics with no lateral expansion is 
used. The second panel is for a uniform jet with sharp 
edges. The two bottom panels are for a Gaussian jet, 
in energy per solid angle, and either a Gaussian or a 
uniform initial Lorentz factor. The viewing angles are 
#obs/#o = 0, 0.5, 1, 1.5, 2, 3, 4, 5, where 80 is the (ini- 
tial) half-opening angle for the uniform jet (two top pan- 
els) and the core ang le (8 C ) for the Gaussian jet (two 
bottom panels). Cfrom lEichler fc GranotJl20o'r?) . 

flat early part of the light curve of XRGRB 041006 
with (0 obs - O ) ~ O.150Q and XRF 030723 with 
(0nhs - O ) ~ 0o l|Granot. Ramirez-Ruiz fc Pernal 
120051 see Figure 12*5)) . as naturally expected in this 
model, since a larger angular displacement outside 
the edge of the jet should result in a softer prompt 
emission. There is a good case for a viewing angle 
slightly outside the edge of a uni form sharp edged jet 
also i n the case of GRB 031203 (|Ramirez-Ruiz et alJ 
I2005L see Figure HO). 

This expectation arises, however, assuming that 
the regions of prominent gamma-ray emission and 
afterglow emission coincide. If this assumption is re- 
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Fig. 25. Tentative fits to the afterglow light curves 
of X-ray rich GRB 041006 (upper panel) and X- 
ray flash XRF 030723 (lower panel) for a uniform 
sharp edged jet with a constant half-opening angle 
8o, viewed from an angle Sobs slightly outside of its 
edge, with (6 ohs - 9 ) ~ O.150 o for XRGRB 041006 
and (fl obs - 6o) ~ 6>o for XRF 0307 23 (for details see 
iGranot. Ramirez-Ruiz fc Pernal l2005t) . In both panels 
the upper curve is optical (R-band) and the lower curve 
is X-rays (0.5-6 keV for XRGRB 041006 and 0.5-8 keV 
for XRF 030723), while the inset shows the optical and 
X-ray spectral slopes. 

laxed l|Eichler fc Granotll2006|) . then an initially flat 
light curve may appear in hard and bright GRBs, 
for lines of sight along which there is bright gamma- 
ray emission but hardly any afterglow emission. In 
this picture viewing angles outside the region within 
the jet with bright afterglow emission can naturally 
account for the early fl at part of the X-ra y after- 
glows detected by Swift (|Nousek et al.ll2006|) . These 
early flat decay stages cannot be attributed to view- 
ing angle effects in the USJ model, where there is 
always significant afterglow emission along the line 




9 0bs =[ 0.5,1 ,1 .25,1 .5,1 .75,2,2.25,2.5]6 
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Fig. 26. A tentative fit to the afterglow light curve of 
GRB 031203 with a uniform sharp edged jet of half- 
opening angle 6o — 5° for di fferent viewing angles ff b s 
from its symmetry axis (from lRamirez-Ruiz et alll2005l) . 

of sight. 15 In scenarios where a break in the light 
curve is caused by the jet structure and/or dynam- 
ics, the break is expected to be largely achromatic. 

3.6. Orphan Afterglows 

For a jet with reasonably sharp edges, the prompt 
GRB emission becomes very dim at viewing angles 
significantly outside the edge of the jet (rA0 o bs ^> 
1), and it can therefore be detected only from 
within or slightly outside the initial jet aperture, 
ro(#obs — do) £ a few. From larger viewing angles 
the prompt emission would not be detected. Dur- 
ing the afterglow, however, the jet decelerates and 
its beaming cone widens with time, so that the af- 
terglow emission at late times can become detectable 
out to much larger viewing angles. Such events, with 
no detected prompt GRB emission but with detected 
afterglow emission in lower frequencies at later times, 
are called orphan afterglows. While no such orphan 
afterglow has been detected to date, their potential 

15 With the possible exception of very large viewing angles, 
$obs > 0max > 0.5 rad, if the jet has a sharp outer edge at an 
angle max , however in this case the event would be very dim, 
and this cannot account for most of the Swift GRBs with an 
early fiat decay phase. 
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for constraining the degree of collimation of the jet 
has been realized early on l)Rhoadslll997h . 

In addition to "off-axis" orphan afterglows, where 
the prompt emission is not detected due to a view- 
ing angle outside of the jet, there can also be "on- 
axis" orphan afterglows, where there is significant 
afterglow emission along the line of sight from rela- 
tively early times, but for some reason the prompt 
gamm a-ray emission along the line of sight is ver y 
weak (jrluang. Dai fc Lul2002tlNakar fc Pirarl2003j) . 
This could occur, for example, if the initial Lorentz 
factor along the line of sight is not large enough to 
avoid excessive pair production (the "compactness" 
problem), T < 100, while it is still sufficiently large 
in order for the deceleration time idee to be early 
enough to enable the detection of the afterglow emis- 
sion (typically i d cc < 1 day for T > 10). 

Many works have analyzed the expected de- 
tection rate of orphan afterglows, as well as 
the constraints on the degree of collimation of 
GRB jets from existing surveys, in th e X-ray 
I Woods fc LoerjlT993|Nakar fc Pi ra iJl200fl>. optical 



Dalai. Griest fc PruitJ l2002t iTotani fc Panaitescii 



iRhoad: 



20021: iNakar. Piran fc Granotf 1200 
.20031: iRau. Greiner k Schwartzl l2006l>. and ra- 
dio JPerna fc Loebl Il998t iLevinson et alJ l2002t 
iGal-Yam et all l200fij) . The constraints derived in 
this way are still not very severe, but are nev- 
ertheless becoming increasingly more interesting. 
Future surveys for orphan afterglows in the X-ray, 
optical and radio can help constrain the degree of 
collimation of GRB jets, as well as the jet structure. 
The detection of orphan GRB afterglows would 
provide an important independent line of evidence 
in favor of jets in GRBs. 

3.7. Some Implications of recent Swift Observations 

The ability of the Swift satellite to rapidly and 
autonomously slew toward GRBs and observe them 
in X-rays, UV, and optical, has dramatically im- 
proved our knowledge of the early afterglow emis- 
sion. In particular it provided excellent coverage of 
the early X-ray afterglow, and enable rapid followup 
observations in the optical and NIR by ground based 
robotic telescopes. In the context of jets, there are 
two main new observations which are most relevant. 

The first is the lack of a clear jet break in the af- 
terglow light curve of most Swift GRBs. Even when 
a break in the light curve does sh ow up, it is often 
chromatic i Panaitescu et al1l2006h . i.e. seen in the 
X-rays but not in the optical, in stark contrast with 
the largely achromatic nature that is expected for 
a jet break. The lack of a clear jet break in many 



Swift GRBs can at least in part be attributed to 
the larger sensitivity of Swift compared to previous 
missions, that causes it to detect dimmer GRBs on 
average, which in turn correspond to wider jets with 
a later jet break time and lower flux at that time, 
making it harder to observe a clear jet break. It is 
not yet clear whether this explanation can fully ac- 
count for paucity of clear jet breaks in Swift GRBs. 
Furthermore, the chromatic breaks that are seen in 
the afterglow light curves of some Swift GRBs defi- 
nitely require a novel explanation. 

The second new Swift observation that bears rel- 
evance for the efficiency of the prompt gamma-ray 
emission, e 7 , and for the kinetic energy, Ej-, of the 
jet during the late phases o f the afterglow for differ- 
ent je t structures ffollowing lGranot. K5nigl fc Piraiil 
l2006h . is the flat decay phase in the early 

Pre- 



Swift studies iPanaitescu & KumaJ 


2002llYost et al l 


2003 


iLlovd-Ronning & Zhandl2004 


) found that the 



isotropic equivalent kinetic energy in the the af- 
terglow shock at late times (typically evaluated at 
t = 10 hr), i?k,iso(10 hr), is comparable to the 
isotropic equivalent energy output in gamma rays, 
£7 7 ,i so , i.e. that typically k = -E 7iiso /-E k ,iso(10 hr) ~ 
1. The gamma-ray efficiency is given by e 7 = 
-^7,iso/ (-^y.iso ~\~ ^k, iso,o), where -Ek,iso,o is the initial 
value of -Ek.iso corresponding to material with a suf- 
ficiently large initial Lorentz factor (ro > 10 2 ) that 
could have contributed to the prompt gamma-ray 
emission. This implies a simple relation, e 7 /(l — 
e 7 ) = Kf, where / = E Klso (10 hx)/E i<isofi can be 
estimated from the early afterglow light curve. 

If the flat decay phase in the early X-ray af- 
terglow obser ved by Swift is interpreted as en- 
ergy injection dNousek et al.l2006HZhang et al. 20061 



IPanaitescu et al.ll2006t iGranot fc Kumarll2006|) this 
typically implies / > 10 and therefore e 7 > 0.9. This 
is a very high efficiency for any reasonable model for 
the prompt emission, and in particular for the pop- 
ular internal shocks model. If the early flat decay 
phase is not due to energy injection, but is instead 
due to an increase with time in the afterglow effi- 
ciency, then / ~ 1 and typically e 7 ~ 0.5. This 
is a more reasonable efficiency, but still rather high 
for internal shocks. Such an increase in the after- 
glow efficiency can occur, e.g., if one or more of the 
following shock micro-physics parameters increases 
with time: the fraction of the internal energy in rel- 
ativistic electrons, e e , or in magnetic fields, eg, or the 
fraction £ e of the electrons that are accelerated to a 
relativistic power-law distribution of energies. If, in 
addition, i?k,iso(10 hr) had been underestimated, e.g. 
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due to the assumption that £ e = 1, then 16 k ~ £ e and 
£ e ~ 0.1 would lead to k~ 0.1 and e 7 ~ 0.1. 

The internal shocks model can reasonably ac- 
commodate gamma-ray efficiencies of e 7 < 0.1, 
which in turn imply n < 0.1. Since the true 
(corrected for beaming) energy output in gamma 
rays, E 1 = fbE Jjiso where f b — (1 - cos 9 n ) 



fln/2 ; is clustered around lO ^ 1 erg (|Frail et all 

l200li iBloom. Frail fc Kulkarnl l2003h . this implies 
£ fc (10 hr) = / 6 £ k , iso (10 hr) = E 7 /k > 10 52 erg 
for a uniform jet. For a structured jet with equal 
energy per decade in 9 (e oc 6* -2 ) in the wings, 
the true energy in the jet is larger by a factor of 
1 + 21n(0 max /0 c ) ~ 10, which implies £^(10 hr) > 
10 53 erg in order to achieve e 7 < 0.1. Such energies 
are comparable (for the UJ model) or even higher 
(for the USJ model) than the estimated kinetic en- 
ergy of the Type Ic supernova (or hypernova) that 
accompanies the GRB. This is very interesting for 
the total energy budget of the explosion. 

4. SUMMARY AND CONCLUSIONS 

Our current understanding of GRB jets has been 
reviewed with special emphasis on the jet dynamics 
(UD and structure (SjSJl. The main conclusions are as 
follows. Semi-analytic models predict a very rapid 
sideways expansion of the jet once its Lorentz factor 
drops below the inverse of its initial half-opening an- 
gle ( H2.1(l . Numerical studies, however, show a very 
modest degree of lateral expansion as long as the jet 
is relativistic. Such numerical studies include both 
an intermediate approach where the hydrodynamic 
variables are integrated over the radial profile of the 
jet, which significantly simplifies the hydrodynamic 
equations f H2.2|) . and full hydrodynamic simulations 

(9229- 

The full hydrodynamic simulations are the most 
reliable of these methods. The fact that the result of 
the intermediate method (described in H2.2| l for an 
initially Gaussian jet agree rather well with hydro- 
dynamic simulations of an initially uniform jet with 
sharp edges, lends credence to its results, which show 
that also for a "structured" jet the lateral expansion 
is very small (even smaller than for an initially uni- 
form or Gaussian jet, because of the smaller gradi- 
ents in the lateral direction) and the distribution of 
energy per solid angle remains very close to its initial 
form as long as the jet is relativistic. 



le lEichler h Wa xman ( 200^) have pointed out a degeneracy 
where the same afterglow observations are obtained under the 
substitution (E,n) — > (_E,n)/£ e and (e e ,eg) — > ^ e (e e ,eg) for 
a value of § e in the range m e /m v < £e < 1, instead of the 
usual assumption of §e = 1. 



The afterglow image is expected to be rather uni- 
form at low frequencies (radio) and more limb bright- 
ened at higher frequencies (optical, UV or X-rays). 
Its morphology and the evolution of its size can help 
probe the jet dynamics and structure, as well as the 
external density profile (§ 12.5(1 . 

The observed jet break in the light curve in a 
uniform jet occurs predominantly due to the lack of 
contribution to the observed flux from outside the 
edges of the jet, once they become visible, and the 
very modest lateral expansion does not play an im- 
portant role f H2.6fl . 

The most popular models for the jet structure are 
the uniform jet (UJ) model, where the jet is uniform 
within some finite half-opening angle and has sharp 
edges, and the universal structured jet (USJ) model, 
where the jet has a narrow core and wings where the 
energy per solid angle drops as a power law (usu- 
ally assumed to be an inverse square) with the angle 
from the jet symmetry axis. There are also other jet 
structures that have been discussed in the literature, 
which include a Gaussian jet, a two component jet 
with a narrow uniform core of initial Lorentz factor 
To > 10 2 surrounded by a wider uniform component 
with Fo ~ 10 — 30, and jets with a cross section in 
the shape of a ring (or "hollow cone") or a fan (see 
Figure [ig. 

There are various approaches for constraining the 
jet structure. Statistical studies of the prompt emis- 
sion are not very conclusive yet, while the observed 
combined distribution of jet break times and red- 
shifts appears to disfavor the USJ model (<JOJ). The 
evolution of the linear polarization of the afterglow 
provides interesting constraints on the jet structure 
and the magnetic field configuration in the emitting 
region f ^3.3[) . but it is difficult to obtain tight con- 
straints on the jet structure without making strong 
assumptions about the magnetic field configuration. 

The shape of the afterglow light curves is an im- 
portant and relatively robust diagnostic tool for con- 
straining the jet structure f H3.4f) . It can practically 
rule out a jet with a cross section of a narrow ring 
or a fan, and it constrains the properties of a two 
component jet. The light curves for viewing angles 
slightly outside the (reasonably sharp) edge of a jet 
would initially rise or at least be very flat before they 
join the decaying light curve for lines of sight within 
the (initial) jet aperture ( H3.5|) . This can naturally 
account for such a flat decay phase observed in the 
best monitored pre- Swift X-ray flash (030723) and X- 
ray rich GRB (041006). It can also explain the early 
X-ray afterglow light curves of many Swift GRBs if 
the regions of prominent gamma-ray and afterglow 
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emission do not coincide ijEichler fc Grano tl l2006|) . 
However, this cannot be attributed to viewing an- 
gle effects in the USJ model, while in the Gaussian 
jet model it suggests that the initial Lorentz factor 
significantly drops outside of the jet core. 

Our understanding of the structure and dynam- 
ics of GRB jets may improve thanks to future ob- 
servations. These include a dense monitoring of the 
afterglow emission starting at early times and over a 
wide range of frequencies (radio, mm, NIR, optical, 
UV, X-ray), polarization measurements with good 
temporal coverage both at early times (much earlier 
than the jet break time) and around (i.e. at least a 
factor of a few before and after) the jet break time, 
and measuring the evolution of the afterglow image 
size. Surveys for orphan afterglows f ij3.6(l are al- 
ready beginning to provide interesting constraints on 
the collimation of GRB jets, and future surveys may 
also constrain the jet structure. The detection of or- 
phan GRB afterglows may also provide an important 
independent line of evidence for jets in GRBs. Re- 
cent Swift observations (323 show a paucity of clear 
achromatic jet breaks in the afterglow light curves 
as well as chromatic breaks which challenge existing 
models and call for new explanations. 

I am grateful to Ehud Nakar, Enrico Ramirez- 
Ruiz, and Arieh Kdnigl for useful comments on the 
manuscript. This research was supported by the US 
Department of Energy under contract number DE- 
AC03-76SF00515. 
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